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SUMMARY 


Teleselsmic  estimations  of  the  yield  of  underground  nuclear 
explosions  are  sometimes  carried  out  by  fitting  the  observed  spectra  and 
amplitudes  of  the  teleselsmic  pulses  to  source  functions  taken  from  near- 
field  data  over  a  range  of  known  events.  The  source  functions  can  be 
parameterized  according  to  yield,  depth,  and  material  type  so  that  given 
some  assumptions  about  depth,  which  may  also  be  available  from  teleselsmic 
data,  and  material,  an  estimate  of  yield  can  be  achieved  provided  that 
properties  of  the  propagation  are  known.  It  has  usually  been  assumed  that 
the  attenuation  of  the  propagating  signal  can  be  described  by  a  Q  function 
which  Implies  that  the  mechanisms  are  linear  In  nature.  The  Q  function  Is 
taken  to  be  a  function  of  location,  especially  depth,  and  It  also  may  be  a 
function  of  frequency.  Establishment  of  the  details  of  Q  is  the  subject 
of  considerable  Investigation  which  Is  largely  based  on  teleseismlc  obser¬ 
vations. 


The  assumption  that  propagation  beyond  the  range  of  a  few  hun¬ 
dred  meters  per  klloton  to  the  1/3  power  Is  In  fact  linear.  Is  the  subject 
of  some  controversy.  It  Is  recognized  from  various  experiments  on  attenu¬ 
ation  of  signals  In  earth  solids  that  the  apparent  Q  tends  to  decrease, 
giving  more  absorption,  when  the  strain  exceeds  10-6 .  Since  the  strains 
at  the  radii  where  the  linear  source  functions  are  taken  to  be  valid  are 
such  that  the  strains  are  considerably  greater  than  10"6 ,  It  Is  suggested 
that  the  mild  nonlinear  effects  which  may  exist  for  a  portion  of  the 
propagation  path  may  significantly  Influence  the  yield  estimates  found, 
assuming  a  linear  Q  function  description. 


In  this  report  we  have  studied  the  remarkably  clean  and  complete 
data  set  taken  from  the  free-fleld  measurements  following  the  nuclear 
event  SALMON.  This  5.3  kt  event  took  place  In  a  rather  uniform  salt  dome 
and  the  data  are  the  best  available  for  the  study  of  properties  of  attenu¬ 
ation  In  the  moderate  strain  region.  We  have  selected  a  set  of  six 
records  from  SALMON  and  attempted  to  establish  whether  or  not  they  are 
consistent  with  the  assumption  of  linear  attenuation  descrlbable  with  a  Q 
function.  It  appears  that  the  data  are  Internally  consistent  with  a  Q  of 
about  10,  which  Is  Independent  of  amplitude  even  though  the  amplitude 
changes  over  more  than  an  order  of  magnitude  over  the  records  used  and  the 
strains  are  typically  about  10-3 .  However,  It  must  be  noted  that  other 
data  from  attenuation  of  oscillation  of  salt  rods  Indicates  that  for  small 
strains  the  Q  Is  about  500.  Thus,  while  the  SALMON  data  are  Internally 
consistent  with  a  constant  Q,  It  Is  difficult  to  see  how  they  can  be  made 
externally  consistent  with  other  data.  Apparently  for  strains  between 
10"4  and  10*6  there  must  be  a  transition  In  behavior  which  Is  clearly 
Indicative  of  nonlinear  behavior.  It  may  be  possible  to  clarify  this 
Issue  by  studying  the  data  from  the  COWBOY  series  of  chemical  explosions 
In  salt  since  they  cover  a  set  of  strains  which  Include  values  down  to 
about  10-5  . 

Inclusion  of  nonlinear  effects  In  attenuation  due  to  physical 
mechanisms  such  as  crack  friction  Is  generally  accomplished  by  providing  a 
description  of  the  microscopic  loss  source  and  Integrating  It  over  a  cycle 
to  find  an  effective  Q.  However,  this  Is  not  very  satisfying  since  Q  Is 
Inherently  a  linear  concept.  In  anticipation  of  the  need  for  a  descrip¬ 
tion  of  nonlinear  attenuation,  we  have  considered  some  constitutive  rela¬ 
tions  which  are  of  such  a  simple  form  that  their  effects  can  be  studied 
analytically.  The  descriptions  of  nonlinear  friction  so  obtained  have 
been  structured  so  as  to  enforce  the  conditions  that  the  effective  Q"1 
Increases  In  proportion  to  the  strain  and  that  the  resulting  waveform 
obeys  the  cube-root  scaling  which  Is  observed  In  a  wide  variety  of  data 


from  explosions.  It  is  found  that  a  simple  form  can  be  found  for  a 
constitutive  relation  which  is  consistent  with  the  experimental  facts, 
including  a  tendency  for  the  waveform  to  change  shape  only  slowly  during 
propagation  over  the  range  seen  In  the  SALMON  data.  It  Is  hoped  that  such 
a  description  will  be  useful  In  characterizing  the  expected  transition  to 
small  strain  behavior. 


SECTION  1 
INTRODUCTION 


It  is  generally  recognized  that  the  near-field  pulse  from  under¬ 
ground  explosions  exhibits  a  nonlinear  behavior,  at  least  out  to  a  radius 
for  which  gross  structual  changes  in  the  rock  appears.  Beyond  this  radius 
(- 100  m/kt1 /3  )  there  may  be  subtler  nonlinear  changes  which  influence 
propagation.  For  the  purposes  of  establishing  a  seismic  source  function 
It  is  common  to  establish  a  range  beyond  which  the  pulse  is  taken  as  line¬ 
ar,  as  Indicated  by  experimental  data  Indicating  a  constant  reduced  dis¬ 
placement  potential.  Typically  this  “elastic  radius"  is  taken  at  -200  or 
300  m/kt1/ 3  and  using  a  model  such  as  that  from  Mueller  and  Murphy1, 
relations  between  teleselsic  observations  and  Implied  source  characteris¬ 
tics  are  deduced.  In  carrying  out  such  calculations  it  is  assumed  that 
the  signal  propagates  linearly  with  some  attenuation  Imposed  in  the  form 
of  Q  function.  If  there  Is  significant  contribution  from  nonlinear  behav¬ 
ior  the  resulting  yield  estimates  may  be  In  error.  In  particular,  there 
Is  evidence?  that  the  pulse  from  explosions  may  have  significant  nonlinear 
features  well  beyond  300  m/kt1/3. 

Of  the  large  body  of  data  available  for  the  propagation  of  seis¬ 
mic  signals  In  the  intermediate  strain  regime  corresponding  to  strains 
varying  from  about  10*3  to  about  10*6  the  data  for  salt  provides  the  best 
information.  If  we  combine  data  from  the  SALMON  nuclear  explosion3,  the 
COWBOY  series4  of  field  chemical  explosions  and  Larson's  data  on  labora¬ 
tory  chemical  explosions5  we  have  available  propagation  data  In  the  rele¬ 
vant  strain  regime  for  a  wide  range  of  yields,  distances  and  characteris¬ 
tic  frequencies;  the  range  of  scaled  distances  (s distance  x  yield*1/3)  is 


a  factor  of  about  300  while  the  range  of  characteristic  frequencies  is 
about  four  orders  of  magnitude.  In  addition  to  the  data  from  explosions 
we  may  add  the  laboratory  data  on  the  absorption  of  the  energy  in  small 
oscillations  of  halite  rods  taken  by  Tittmann6. 

When  viewed  in  the  large  looking  at  the  general  trends  of  the 
data  certain  characteristics  appear.  In  this  section  we  shall  describe 
these  characteristics  and  give  the  elements  of  the  data  which  support 
them.  In  Section  2  we  give  the  results  of  a  detailed  analysis  of  data 
from  event  SALMON— the  event  for  which  the  data  best  permits  a  data  set 
analysis  of  the  propagation  of  the  seismic  wave.  We  found  in  distinct 
contradiction  to  our  prior  expectations  that  the  propagation  is  very  well 
described  by  linear  propagation  In  an  inelastic  medium.  This  observation 
is  the  most  important  result  In  the  report.  In  Section  3  we  give  a  non¬ 
linear  constitutive  relation  which  provides  for  propagation  with  proper¬ 
ties  thought  characteristic  of  the  overall  data.  The  detailed  studies  of 
the  propagation  of  the  SALMON  shock  show  that  some  features  which  appear 
to  be  characteristic  of  the  entire  data  set  are  not  characteristic  of  sub¬ 
sets  of  it;  in  particular  the  constitutive  relation  we  shall  give  is  not 
appropriate  for  calculating  the  propagation  of  the  SALMON  shock;  whether 
it  is  useful  for  calculating  the  propagation  of  waves  as  observed  in  other 
parts  of  the  data  set  is  a  question  which  answer  will  require  further 
examination  of  the  data. 

PROPERTIES  OF  THE  DATA 

These  properties  appear  to  be  characteristic  of  the  propagation 
of  waves  in  the  Intermediate  strain  region: 

1)  The  absorption  of  energy  increases  with  ampli¬ 
tude;  one  relation  which  has  been  proposed  to 
express  this  effect  is: 


where  e  Is  the  peak  strain;  the  quantity  Q  will 
be  discussed  below. 

2)  The  propagation  satisfies  cube-root  scaling2; 

that  Is,  the  propagation  Is  Invariant  under  the 
scale  transformations: 

t-  y1/3  ts 
x  ♦  y1  /3  xs 
u  +  y1  /3  us 

where  t$,  x$  and  u$  are  respectively  the  scaled 
time,  scaled  distance  and  scaled  displacement 
while  y  Is  the  event  yield. 

3)  The  propagation,  although  presumably  nonlinear 

(property  1),  nevertheless  provides  a  waveform 
which  changes  only  slowly  with  distance. 

Let  us  first  consider  the  justification  for  item  1.  As  indica¬ 
ted  In  the  next  section,  the  SALMON  shock,  with  measured  strains  in  the 
range  4  •  10"3  -  3  •  10”'*,  propagated  very  closely  like  a  wave  In  an  anelas- 
tic  solid  with  a  Q  of  -  10-20;  laboratory  measurements  by  Tlttmann6  give  a 

value  of  Q  of  -500  for  strains  greater  than  about  10“®  with  the  apparent  Q 

beginning  to  Increase  with  amplitude  when  the  strain  becomes  greater  than 
10-6  .  Tlttmann's  data  were  taken  at  a  frequency  of  the  order  of  500  Hz, 
which  Is  only  somewhat  higher  than  the  frequency  at  which  the  SALMON  data 
appears  to  be  very  good;  thus  the  factor  of  -50  difference  in  Q  between 
the  SALMON  data  and  Tlttman's  data  does  not  appear  to  be  a  frequency 
effect;  presumably  It  is  an  amplitude  effect.  This  argument  provides  the 
strongest  reason  for  believing  In  Item  1. 


In  addition  to  the  argument  just  presented  various  workers  have 
attempted  to  formulate  arguments  In  favor  of  Item  1  by  using  peak  veloc¬ 
ity  data.  Minster  and  Archambeau7  proposed  a  relationship  of  the  form 
given  In  Item  1,  pointing  out  that  such  a  relationship  Is  consistent  with 
the  general  trend  of  the  SALMON-COWBOY  data  and  also  with  Tlttmann's 
laboratory  data.  Later,  the  present  authors8  reached  the  same  conclusions 
based  on  the  same  considerations  but  concluded  that  Q-1  proportional  to 
✓e  or  e2  would  work  about  as  well  as  the  linear  relation  given  In  item  1. 
(See  Appendix  A,  where  this  paper  Is  reproduced.)  The  fact  Is  such  argu¬ 
ments  are  very  hard  to  make  In  any  convincing  way.  The  problem  Is  that 
the  decrease  of  peak  amplitude  depends  on  a  combination  of  absorbed  energy 
and  charged  waveform  and  the  freedom  this  mix  allows  In  possible  propaga¬ 
tion  laws  causes  the  peak  amplitude  data  to  be  only  a  weak  constraint  on 
the  underlying  physics.  In  spite  of  these  difficulties  and  reservations 
we  shall,  for  the  moment,  adopt  the  attitude  that  item  1  Is  a  characteris¬ 
tic  of  the  propagation  and  shall  attempt  to  formulate  a  constitutive  rela¬ 
tion  which  leads  to  solutions  consistent  with  It. 

The  major  data  supporting  Item  2  Is  presented  In  Figure  1.  In 
the  figure  the  peak  particle  velocity  is  plotted  against  scaled  range 
(xy-i/3)  f0r  a  number  of  laboratory  experiments,  field  experiments  with 
chemical  explosives  and  the  nuclear  explosion  SALMON.  It  appears  that  we 
can  write 


v(yx)  »  V(x,y-1/3 )  »  v(x$)  (1) 

In  this  relation  we  have  suppressed  the  time  dependence;  the  fact  that.  In 
the  strain  regime  we  are  concerned  with  the  signals  travel  with  the  sound 
speed  of  the  medium  Independent  of  the  yield  implies: 


v(y,x,t)  ■  Vfxy-1  /3  ,ty-1  fi  )  •  v(x  .t  ) 


(2) 


Scaled  range,  m/kt,/3 


Figure  1.  Peak  particle  velocities  from  explosions  in  salt  (Reference  5) 


Integrating  equations  (2)  with  respect  to  time  we  find 


u(y,x,t)  »y 1/3  /  7(xs,ts)  dt$  ■  yI/3u$  (3) 

where  u  Is  the  particle  displacement.  Thus  we  obtain  the  scaling  laws  of 
Item  2. 


The  argument  just  presented  Is  again  based  on  peak  velocity  data 
(Figure  1).  In  this  case  however  we  believe  that  the  argument  Is  much 
more  persuasive  than  for  Item  1.  The  various  experiments  represented  In 
Figure  1  corresponded  to  Initial  conditions  which  differed  to  at  least 
some  degree  while  the  scaling  rule  appears  to  be  very  closely  satisfied. 
It  seems  extremely  likely  that  the  scaling  rules  of  item  2  are  a  property 
of  the  underlying  dynamical  equations  rather  than  a  combination  of  dynam¬ 
ics  which  does  not  satisfy  scaling  and  particular  Initial  conditions.  We 
shall  assume  that  our  constitutive  relations  should  lead  to  dynamical 
equations  which  satisfy  scaling. 

The  scaling  property  of  Item  2  Is  often  referred  to  as  "fre¬ 
quency  Independence  of  Q";  let  us  understand  why  that  Is:  If  we  Fourier 
transform  the  velocity  we  find: 

v(y,x,t)  *  /  e1“tv(y ,x,t)dt 

-  yl/3  Jei<0stV(xs.tx)dts  (4) 


y1/3  v  (x.w.) 
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Let  us  now  assume  that  the  velocity  propagates  from  one  point  to  the  next 
according  to: 


vfy.Xj  *,) 


v(y,X!  ^») 
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which  Is: 
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2Qc 
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(6) 


We  assume  that  the  product  Qc  Is  not  a  function  of  distance  or  yield; 
therefore,  it  cannot  be  a  function  of  u  either,  else,  the  R.H.S.  would  not 
be  a  function  of  scaled  variables  whereas  the  L.H.S.  Is.  Observational ly 
c  Is  nearly  Independent  of  frequency;  therefore  Q  must  be  also. 

We  thus  see  that  for  a  linear  theory  scaling  Implies  that  Q  is 
independent  of  frequency.  For  a  nonlinear  theory  there  are  two  problems 
with  this  statement.  The  first  Is  that  for  general  nonlinear  propagation 
there  is  no  obvious  definition  of  Q.  The  second  is  that  for  special  types 
of  nonlinear  propagation  where  there  Is  an  obvious  definition  of  Q  it  may 
no  longer  be  true  that  scaling  Implies  that  Q  Is  frequency  independent,  as 
we  shall  see  in  an  example  given  In  Section  3.  In  thinking  about  non¬ 
linear  propagation  we  must  remember  that  the  data  of  Figure  1  implies 
scaling,  not  frequency  Independence  of  Q. 

Item  3  Is  meant  to  be  only  a  qualitative  statement,  the  degree 
to  which  It  Is  true  can  be  judged  by  comparing  Figures  2a  and  2b  which 
show  the  shock  from  the  same  event  at  two  ranges  for  a  representative 
case. 
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SECTION  2 

SALMON  DATA  AND  NONLINEAR  BEHAVIOR 


In  this  section  Me  will  examine  some  data  from  SALMON,  which  are 
among  the  cleanest  available.  In  order  to  see  If  there  Is  not  definitive 
evidence  for  nonlinear  features  beyond  the  "elastic  radius."  This  will  be 
done  by  examination  of  a  series  a  velocity  records  from  subsurface  points 
to  see  If  they  are  consistent  with  strictly  linear  behavior  since  the 
range  of  linear  effects  Is  relatively  limited.  We  shall  attempt  to  deter¬ 
mine  If  the  records  exhibit  features  which  cannot  be  understood  In  terms 
of  linear  attenuation  with  an  absorption,  which  Is  possibly  a  function  of 
frequency.  Since  this  Is  a  form  to  which  any  linear  mechanism  can  be 
reduced.  It  will  serve  to  bound  the  range  of  linear  behavior.  Basically 
If  the  data  are  consistent  with  attenuation,  which  Is  Independent  of  amp¬ 
litude  but  possibly  a  function  of  frequency.  It  may  allow  the  possibility 
that  the  propagation  mechanism  Is  linear  without  ruling  out  a  nonlinear 
effect.  However  It  can  serve  to  rule  out  a  wide  range  of  possible  non¬ 
linear  mechanisms.  In  fact  this  Is  exactly  the  outcome  of  this  section. 

NEAR-FIELD  DATA 

While  there  exists  near-field  data  from  shots  In  a  variety  of 
media?,  by  far  the  most  consistent  set  are  from  explosions  In  salt2*3. 

Salt  Is  a  relatively  uniform  medium  for  which  a  variety  of  experiments 
ranging  from  SALMON,  through  the  COWBOY  series?*  to  small  laboratory  explo¬ 
sions  as  reported  by  Larson5.  These  data  provide  the  basis  for  many  atti¬ 
tudes  held  concerning  properties  of  near-field  data  so  they  will  be 
reviewed  at  this  time.  The  experiments  of  Larson  for  small  chemical 
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explosions  In  pressed  salt  have  provided  pulses  over  a  scaled  range  from 
10  m/kt1/3  to  200  m/kt1/3.  The  range  of  frequencies  covered  was  from 
about  104  to  106  Hz  and  the  ratio  of  peak  particle  velocities  to  compres- 
slonal  sound  speed  (which  Is  comparable  with  the  strain)  went  from  about 
10" 1  to  less  than  10" 3,  which  by  most  standards  would  suggest  that  the 
response  was  nonlinear.  Yet,  by  performing  a  direct  superposition  experi¬ 
ment  with  a  pair  of  simultaneous  explosions.  It  was  found  that  the  result¬ 
ing  response  was  consistent  with  direct  addition  of  the  two  pulses  as 
would  be  expected  from  a  linear  medium.  Still  It  Is  not  clear  just  how 
nonlinear  effects  would  be  manifest  In  this  experiment  without  knowing  the 
character  of  any  nonlinear  behavior.  That  Is,  the  apparent  success  of 
superposition  for  pulses  with  large  strains  may  not  directly  negate  the 
possibility  of  any  sort  of  nonlinear  behavior.  The  Cowboy  series  of 
chemical  explosions  had  a  range  of  yields  from  10  to  2000  pounds  of  TNT, 
some  of  which  were  carried  out  In  cavities  for  decoupling  tests.  The 
scaled  ranges  of  these  experiments  were  from  200  m/kt1/3  to  3000  m/kt1/3 
and  the  corresponding  peak  strains  were  from  a  few  times  10-4  to  about 
10"5.  The  SALMON  event,  «t  5.3  kt,  took  place  In  a  natural  salt  dome  and 
a  comprehensive  set  of  measurements  were  taken,  both  at  surface  and  sub¬ 
surface  locations.  These  subsurface  measurements  Included  scaled  ranges 
from  100  m/kt1/3  to  about  425  m/kt1/3,  which  provided  peak  strains  from 
about  4. 10" 3  to  about  3 »10" 4  with  a  significant  frequency  content  from 
about  1  Hz  to  more  than  102  Hz.  The  SALMON  data  are  well  known  for  their 
Internal  consistency  and  their  correspondence  with  the  other  salt  data. 
Peak  velocity  data  from  the  salt  shots  were  shown  In  Figure  1  as  a  func¬ 
tion  of  scaled  range.  The  remarkable  result  Is  that  the  data  from  a  huge 
range  of  yields  tends  to  nearly  fall  on  a  straight  line  Indicating  a  power 
law  behavior  with  an  exponent  of  about  -1.9.  This  contrasts  with  a  value 
of  -1  which  would  be  expected  for  elastic  behavior.  Furthermore,  It  Is 
remarkable  that  the  simple  scaling  with  yield1/3  seems  to  hold  over  the 
great  range  of  strains.  As  Trullo2  has  noted,  the  nuclear  event  appears 
to  scale  slightly  differently  from  chemical  explosions  but  this  Is  expec¬ 
ted  based  on  the  different  time  scale  of  their  source  functions  and  can  be 
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accounted  for  by  simply  altering  yields  by  a  multiplicative  constant  to 
convert  to  chemical  yield.  It  should  be  noted  that  the  range  of  yields 
over  which  simple  scaling  holds  Includes  strains  which  are  expected  to 
give  nonlinear  behavior.  However,  It  has  not  been  cleanly  demonstrated 
that  the  deviation  from  r-1  behavior  Is  necessarily  a  nonlinear  effect. 
In  particular.  It  has  been  observed  that  linear  but  Inelastic  (l.e., 
anelastlc)  behavior  may  provide  such  results.  Trull o  has  noted  that  the 
SALMON  data  for  decay  of  peak  velocity  with  range  are  consistent  with  an 
effecttve  Q  of  about  5.  Here  Q  Is  found  by  considering  and  amplitude 
decay  with  range.  In  addition  to  geometrical  divergence,  of  the  form 

exp(-«r/2cQ) 

Note  that  simple  scaling.  In  conjunction  with  linearity.  Indicates  that  Q 
must  be  Independent  of  frequency  and  Independent  of  amplitude.  To  some 
degree  both  of  these  suggestions  are  subject  to  experimental  examination. 
The  work  of  Larson  on  explosions  In  laboratory  pressed  salt  Indicate  that 
while  superposition  appears  to  hold,  values  of  Q  of  12  and  25  with 
Increasing  range  suggest  that  a  constant  Q  cannot  describe  the  results. 
Furthermore,  laboratory  experiments  by  Tlttmann  on  decay  of  cyclic  motion 
Induced  In  salt  samples  Indicate  that  for  strains  below  10-6  a  value  of  Q 
of  -500  Is  appropriate  but  that  for  larger  strains  the  value  of  Q  (more 
properly,  the  effective  Q)  decreases  Indicating  a  nonlinear  attentuatlon. 
While  these  various  experiments  took  place  In  different  frequency  regimes, 
there  Is  no  known  frequency  dependence  which  can  tie  them  together  under  a 
linear  theory.  The  SALMON  data  provide  a  remarkably  consistent  set  which 
may  allow  for  drawing  clean  conclusions  concerning  the  character  of  near- 
field  propagation  without  the  necessity  for  Introduction  of  any  assump¬ 
tions  about  the  nature  of  scaling  among  unrelated  events.  Perret3  has 
presented  these  data  In  comprehensive  form  and  provided  fits  without 
regard  to  any  particular  theoretical  assumptions.  These  data  were  further 
used  by  Trullo  In  his  analyses  of  the  degree  to  which  they  Indicate  scal¬ 
ing.  In  this  report  we  will  consider  the  extent  to  which  the  SALMON  data 
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can  shed  light  on  the  Issue  of  linear  versus  nonlinear  behavior  by  explor 
Ing  the  consequences  of  the  assumption  of  linear  behavior. 


SALMON  DATA  REDUCTION 

Perret  has  provided  plots  of  the  SALMON  near-field  pulses  for 
accelerations  as  well  as  velocities  and  displacements.  The  original  data, 
as  digitized  at  Sandla,  are  available  and  they  were  kindly  supplied  to  us 
by  J.  Trullo  and  N.  Perl  of  Applied  Theory,  Inc.  The  use  of  the  original 
data  Is  preferred  to  tracing  Perret 's  curves  since  Perret  has  made  certain 
assumptions  In  correcting  the  data  to  force  zero  acceleration  and  velocity 
at  late  time.  In  his  report,  Perret  did  not  Indicate  the  sorts  of  correc¬ 
tions  which  are  required  and  Inasmuch  as  the  character  of  the  corrections 
may  well  Influence  the  outcome  of  spectral  methods  we  felt  It  best  to 
return  to  the  original  data  In  order  to  have  control  over  all  operations 
on  the  data.  The  general  nature  of  the  nonphysical  character  of  the  raw 
data  Is  Indicated  by  Figure  3  which  shows  a  once  Integrated  acceleration 
record  to  give  velocity. 

It  will  be  noted  that  two  problems  are  evident.  First,  the 
velocity  not  only  does  not  go  to  a  constant  but  It  Increases  at  a  constant 
rate.  This  Is  as  result  of  a  post-shot  non-zero  baseline  for  the  accele¬ 
ration  Instrument.  Evidently  the  pulse  from  the  explosion  disturbed  the 
Instrument  zero.  This  can  be  corrected  for  by  altering  the  baseline  at 
late  times.  There  Is  a  problem  beyond  this  because  even  when  the  late 
time  acceleration  Is  taken  as  zero,  the  late  time  velocity,  while  cons¬ 
tant,  Is  not  zero.  This  Is  generally  thought  to  be  a  result  of  clipping 
of  the  acceleration  peak  due  to  Inadequate  Instrument  response  or  to 
Inadequate  bandwidth  In  recording.  It  Is  Impossible  to  correct  the  dif¬ 
ficulties  In  a  unique  manner.  However,  It  appears  to  be  possible  to  do  so 
In  such  way  that  the  effect  on  conclusions  drawn  from  the  result  will  not 
be  Important.  We  have  chosen  to  correct  the  data  by  the  following  proce¬ 
dure  as  illustrated  In  Figure  4. 
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Figure  3.  Raw  velocity  record  f row  SALMON  (E14C27AR)  with  straight  line  used 
to  correct  for  baseline  shifts. 
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First  the  raw  velocity  records,  which  are  given  In  the  appendix,  were  fit 
at  late  time  with  a  straight  line  from  the  latest  time  shown  (1.6  sec) 
back  to  the  time  of  the  peak  In  such  a  way  as  to  give  a  behavior  past  the 
peak  which  was  consistent  with  the  results  shown  by  Perret.  The  velocity 
before  the  peak  was  then  altered  by  multiplication  of  each  data  value  by 
the  ratio  of  the  new  peak  value  to  the  old  peak  value.  This  assured  that 
the  velocity  would  be  continuous  although  with  a  mild  discontinuity  in  the 
first  derivative  exactly  at  the  peak.  It  is  Impossible  to  detect  this 
discontinuity  by  merely  looking  at  the  resulting  plotted  velocities.  As  a 
result  we  expect  that  the  effects  of  this  correction  will  be  seen  primar¬ 
ily  In  the  high  frequency  end  of  the  spectra  and  generally  this  will  not 
be  the  region  of  Interest.  In  addition  to  altering  the  raw  data  records 
for  obvious  Instrumentation  problems,  we  have  at  the  same  time  Introduced 
a  geometrical  factor  which  provides  an  estimate  of  the  spherical  radial 
velocity  from  each  record.  Since  all  six  of  the  records  chosen  for  mani¬ 
pulation  are  of  horizontal  radial  motion,  it  Is  a  straightforward  matter 
to  assume  that  the  motion  was  strictly  along  a  spherical  radius  and  pro¬ 
duce  the  spherical  motion  as 

v(R)*v(r)/cos(e  ) 

where  0  Is  the  angle  with  respect  to  horizontal  of  the  line  through  the 
Instrument  from  the  explosion.  In  several  cases  there  do  exist  records  of 
the  vertical  component  of  the  motion  at  the  same  station  and,  in  prin¬ 
ciple,  one  could  either  combine  these  with  the  radial  component  in  a 
vector  fashion  or  use  them  directly  to  estimate  the  spherical  component. 
However,  examination  of  a  few  cases  Indicated  that  nothing  was  gained  in 
terms  of  consistency  by  using  such  methods  so  we  have  worked  exclusively 
with  the  radial  components. 


The  six  records  selected  for  study  include  those  at  ranges  from 
166  to  660  meters.  The  records  were  selected  so  as  to  have  as  large  a 


set  of  ranges  as  possible  while  also  having  them  be  as  internally  consis¬ 
tent  as  possible.  The  full  set  of  original  data  records  to  which  we  had 
access  is  shown  in  Appendix  B,  along  with  indications  of  the  locations  of 
the  Instruments  on  which  they  were  recorded.  Consistency  was  established 
by  taking  records  whose  peak  velocities  fell  along  a  smooth  curve  when 
plotted  against  range.  While  there  is  no  real  reason  to  believe  that  any 
one  of  the  records  is  more  desirable  than  some  which  were  rejected,  our 
use  of  spectral  ratio  methods  require  us  to  avoid  examples  which  will 
obviously  lead  to  unreasonable  results  when  records  are  used  in  pairs. 
The  records  selected  are  indicated  in  Table  1  along  with  some  of  their 
properties. 


Table  1.  SALMON  records  used. 


Record 

Range 

(m) 

VVr 

VR 

peak 

(m/s) 

mi 

E-14C-27-AR 

166 

1.0 

13.8 

14.0 

E-14C-22-UR 

225 

1.36 

8.0 

8.2 

E-14-20-AR 

276 

1.61 

5.1 

5.5 

E-6-27-AR 

318 

1.0 

3.75 

3.9 

E-14C-39-AR 

402 

2.40 

2.5 

2.7 

E-11-34-AR 

660 

1.06 

1.1 

0.99 

*Perret's  values 


The  original  records  extend  out  to  about  five  seconds  past  the  explosion. 
It  was  found  that  after  about  1.5  seconds  the  records  showed  no  further 
significant  contribution  so  we  truncated  the  data  at  1.6  seconds  giving  a 
convenient  8192  data  points  for  each  record.  The  corrected  velocity 
records,  some  of  which  were  taken  from  acceleration  data,  are  shown  on 
Figure  5. 


When  plotted  on  the  same  scale  it  becomes  difficult  to  fully 
appreciate  the  features  of  the  more  distant  examples  so  the  same  data  are 
plotted  again  in  Figure  6,  suitably  normalized  to  a  unit  peak  amplitude. 

It  appears  from  the  results  in  Figure  5  that  the  pulse  progres¬ 
ses  smoothly  to  large  radii  with  a  decrease  in  amplitude  which,  it  will  be 
shown.  Is  approximately  like  radius  to  the  -1.9  power.  At  the  same  time 
Figure  6  indicates  that  the  pulse  shape  remains  fairly  stable  but  it  shows 
a  mild  tendency  to  smooth  out  the  sharp  peak  and  to  broaden  the  width. 
Both  of  these  features  suggest  that  the  higher  frequency  components  are 
being  more  rapidly  attenuated  than  the  low  as  would  be  appropriate,  for 
example,  with  a  constant  Q  should  such  a  description  be  suitable.  Beyond 
this  we  notice  an  Inital  ramp  on  each  pulse  which  strengthens  relative  to 
the  peak  amplitude  as  the  radius  Increases.  This  feature  is  generally 
attributed  to  shear  failure,  which  occurs  when  the  strain  becomes  suffic¬ 
iently  large,  providing  a  reduced  compressional  velocity.  Note  that  the 
leading  edge  of  the  pulses  propagate  at  about  4.7  km/s  while  the  peak 
propagates  at  about  3.7  km/s.  In  the  following  discussions  we  shall  not 
make  further  reference  to  this  nonlinear  feature. 

In  order  to  establish  the  possible  nonlinear  character  of  the 
SALMON  data  we  shall  adopt  the  position  of  exploring  the  consequences  of 
the  assumption  of  anelastlc  behavior  as  expressed  in  terms  of  a  Q  function 
which  we  will  attempt  to  evaluate.  If  this  attempted  description  leads  to 
contradictions,  such  as  an  effective  Q  which  depends  upon  amplitude,  we 
will  attribute  the  effects  to  nonlinear  behavior.  Basically  we  shall  use 
the  data  to  make  estimates  of  Q  over  a  range  of  frequencies  and  ampli¬ 
tudes.  Since  we  are  dealing  with  data  exclusively  from  a  single  event,  no 
assuptlons  about  the  scaling  behavior  are  necessary  and  no  such  informa¬ 
tion  will  be  gained. 


SCALED  VELOCITY 


Figure  6.  SALMON  records  from  Figure  5  as  plotted  with  common  peak  value  to  <-*— 

Illustrate  the  waveforms. 


A  strict  Interpretation  of  anelastlc  behavior  would  require  that 
we  find  a  function  Q,  dependent  generally  upon  frequency  as  well  as  a 
phase  velocity  c  which  Is  consistent  with  Q,  so  as  to  enforce  the  require¬ 
ment  of  causality.11  However,  moderate  values  of  Q  are  consistent  with  a 
nearly  constant  c,  at  least  over  the  range  of  frequencies  which  are  avail¬ 
able  to  experimental  verification.  This  feature  can  also  be  noted  from 
experimental  data  by  looking  for  dispersion  In  the  implied  p-wave  phase 
velocities.  The  SALMON  data  cannot  be  made  to  establish  the  frequency 
dependence  of  c  due  do  substantial  scatter  In  the  phase  of  the  Fourier 
components  of  the  signal. 

The  basic  structure  of  the  equations  of  motion  for  a  spherical 
pulse  propagating  In  an  Isotropic  medium  can  be  expressed  as 


p  *1“  »  L.(i_(r2G  J)  -  i, 

r2  «p  ^  r,J  r 


3t2  ^2  a  r  I*  08 

where  the  stress  components,  or,  are  related  to  the  strain  components, 
which  are  provided  by  the  radial  displacement  u  through 


Of  course  It  Is  necessary  to  have  or  find  some  constitutive  relation 
between  e  and  a  In  order  to  fix  the  propagation  of  a  pulse.  If  a  linear 
relation  between  e  and  a  is  taken,  a  general  form  which  is  convenient  to 
work  with  Is  as  follows 

opr  *  (x+2u)  /  m(t-r  )crr(x  )  dr 
+  2x  f  m(t-r  )e00  (x  )  dr 
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+  X  f  m(t-r  )e  rr(r  )  dt 


a 
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where  ra(t)  Is  the  derivative  of  the  relaxation  function.  If  m(t)  Is  a 
delta  function  the  relations  reduce  to  the  standard  elastic  result  with 
the  Lame  constants  x  and  y  .  If  we  Introduce  the  definition  of  the 
reduced  displacement  potential,  as 

u(r,t)  ■  (r,t)/r)  (9) 

3r 

the  equation  of  motion  can  be  written  In  compact  form  In  the  frequency 
domain  as 


1  !_(i  u2*  ♦  ;<•)(»♦  2u>  iL;>)  -  o  do) 

3  r  r  g  f2 

where*  will  generally  be  a  function  of  both  r  andw.  From  this  form  we 
>  see  that  the  advancement  of  4  can  be  expressed  as 

*(r,t)  *  L_  /  *(0,w)  e"1ut  e1“r/°^do  (11) 

2>r 

where  the  apparent  p-wave  velocity  Is  given  by 

a2  -  (X+2y)/p  (12) 

and  *(0,u)  Is  the  Fourier  transform  of  the  *  at  r*0. 

Kjartannson12  has  shown  that  there  Is  a  choice  for  m(t)  which 
will  give  a  Q  which  Is  exactly  constant  and  which  provides  a  phase 
velocity 


which  Is  nearly  constant  over  a  wide  range  of  frequency  provided  Q  Is  not 
too  small.  This  choice  for  m  Is 


m(w  )  -  |2Y  exp(1rrsgn(u))  (13) 

'«o 

where  y  and  uo  are  arbitrary  real  parameters.  Thus  a/ m  Is  the  complex 
p-wave  velocity.  The  observed  pulse  propagation  velocity  can  be  used  to 
fix  the  value  of  ai0  and  the  parameter  y  can  be  expressed  In  terms  of  a 
constant  Q  by  noting  that  the  phase  relation  between  e  and  s  gives  Q  which 
can  then  written  as 

(T1  *tan(mr)  (14) 


so  we  can  write 


«*(r,t)  «i_/  du  ♦(Ow-)  exp(-1w(t-£))  •  exp(-u»r  sgn(u )/(2cQ)) 

at  C 

(15) 

where 

c  *  Colw-I7/  cosUy/2)  (16) 

'uo  1 

so  than  the  displacement  will  be  advanced  by 


u(r,t)  «  -1-  /  da  u(a^j)  exp( -1u(  t-  -t=i)  -u(r-a)  sgn(u  )/(2cQ)) 


Jol  -  d>  sqn(m)  _  i_ 

x  cr  2cQr  r2 

JkiL  -  in  sqn(iri)  _  L_ 

ca  2cQa  £ 


(17) 


While  Equation  (15)  Is  specific  to  the  Kjartannson  assumptions. 
It  provides  a  generic  form  which  Is  often  taken  to  be  an  Indication  of 


the  attenuation  which  Is  associated  with  Q.  That  Is,  we  see  that  Q  forces 
an  exponential  decay  with  range  of  the  RDP  by  a  factor 


exp(  -  «-C_) 
*  2cQJ 
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and  this  provides  a  common  definition  of  Q  which  Is  applied  to  experimen¬ 
tal  data  by  taking  the  phase  velocity  c  to  be  constant,  normally  at  the 
nominal  speed  of  the  pulse.  Again,  this  Is  possible  only  because  c  Is  at 
most  a  slowly  changing  function  of  frequency  so  that  the  pulse  retains  a 
form  which  does  not  change  greatly  with  range.  The  data  from  SALMON  as 
shown  In  the  Figure  2  do  retain  the  pulse  shape  to  such  a  degree  that  an 
Initial  presumption  of  constant  c  can  be  made.  In  fact  In  the  following 
analysis  we  will  estimate  the  product  of  cQ  and  then  provide  Q  subject  to 
particular  values  of  c. 

Manipulating  the  relations  among  the  ROP,  Q  and  the  velocity  In 
the  frequency  domain  provides  an  operational  definition  of  the  experimen¬ 
tal  value  of  Q  through  the  equation 


la  _  u _ _  1 

pxl^.(l>(r*a))  .  I v(f VM ) j  .ca  2cQa 

2cQ  v(a^ii)  Jsa. 

cr  2c?F  ,4 


(18) 


which  provides  an  Implicit  expression  for  Q  from  the  data.  In  practice  It 
Is  found  that  the  second  term  on  the  right  hand  side  of  (18)  Is  nearly 
Independent  of  Q  for  Q  not  too  small  so  that  one  can  write 


This  Is  the  usual  “definition"  of  Q  from  experimental  data  although  It 
does  not  exactly  conform  to  the  proper  definition13.  Note  that  at  large 
ranges  the  1/r2  term  Is  negligible  so  the  expression  becomes  yet  simpler. 
However  In  the  following  we  will  take  Equation  (19)  as  the  definition  of  Q 
from  the  data.  One  could  solve  more  accurately  using  (18)  but  our  calcu¬ 
lations  Indicate  that  there  Is  no  significant  difference  between  (18)  and 
(19)  for  the  SALMON  case.  Based  on  the  motion  of  the  peak  of  the  SALMON 
pulse  motion  we  take  c*3.7  km/s  as  the  constant  sound  speed. 

Using  the  six  data  records  given  in  Figure  5,  one  can  form  a  set 
of  five  sequential  estimates  for  Q  at  ever  increasing  ranges.  The  Q  thus 
obtained  Is  explicitly  a  function  of  frequency  and  based  on  a  minimum  sam¬ 
ple  spacing  of  200  microseconds  over  8192  samples,  the  frequency  range  so 
obtained  Is  from  .6  to  2500  Hz.  However,  we  find  that  the  high  frequency 
end  of  the  scale  Is  not  well  defined  by  the  data,  due  In  part  to  the  mani¬ 
pulations  to  avoid  Improper  late-time  behavior  and  in  part  to  noise  in  the 
data.  As  a  result  we  present  only  that  portion  of  the  frequency  regime 
below  500  Hz  with  the  suggestion  that  perhaps  only  the  region  below  100  or 
20Q  Hz  is  meaningful.  Figures  7  through  11  show  the  estimated  Q(<o )  found 
from  successive  records  at  166,  225,  267,  318,  402,  and  660  meters  from 
the  explosion.  Figure  12  gives  an  average  over  the  entire  range  by  using 
only  the  closest  and  most  distant  records.  It  appears  that  the  record  at 
318  meters  Is  not  very  consistent  with  the  others  but  when  the  results  in 
Figures  9  and  10  are  discarded  some  features  of  the  Q  estimates  become 
clear.  First  comparing  Figures  7  and  11  should  Indicate  the  extent  of 
amplitude  dependence  of  this  effective  Q  and  so  point  to  nonlinear  beha¬ 
vior.  However  these  two  estimates  of  Q  are  remarkably  similar  and  they 
also  are  very  much  like  the  average  estimate  for  0  from  Figure  12,  which 
also  gives  the  RVP  at  the  extreme  radii. 

In  all  cases  the  Q  estimates  at  frequencies  above  about  200  Hz 
show  such  great  scatter  that  they  are  probably  meaningless.  However  below 
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Figure  9.  Estimate  of  Q(«)  from  SALMON  records  at  276  and  318  meters. 
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Figure  10.  Estimate  of  Q(«)  from  SALMON  records  at  318  and  402  Meters 


Figure  12a.  Estimate  of  Q(M)  from  SALMON  records  at  166  and  660  meters. 
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200  Hz  the  results  are  remarkably  consistent  in  that  they  show  a  fairly 
constant  Q  down  to  a  few  tens  of  Hz,  followed  by  a  smooth  rolloff  to  zero 
at  low  frequency.  One  should  keep  in  mind  that  the  pulses  from  SALMON 
have  a  typical  period  of  perhaps  0.2  seconds  with  structure,  which  appears 
significant  as  small  as  0.01  seconds.  Thus  the  portion  of  the  spectrum 
which  is  highly  visibly  ranges  from  less  than  5  to  perhaps  100  Hz.  It  is 
thus  of  interest  to  see  if  the  apparent  frequency  dependence  of  Q  stems 
from  features  of  the  pulses  which  are  evident  in  the  time  domain.  We  have 
checked  this  by  reconstructing  the  outer  five  pulses  from  the  inner  one 
making  two  different  assumptions  about  the  behavior  of  Q.  First  consider 
the  case  of  constant  Q  with  c  allowed  to  change  with  frequency  according 
to  the  results  from  Kjartannson  shown  before.  For  this  purpose  we  take  a 
value  of  Q=10  and  c  at  1  Hz  of  3.7  km/s.  Now  the  result  in  Equation  (15), 
when  applied  to  velocity,  can  be  used  to  propagate  the  observed  at  166 
meters  subject  to  exactly  constant  Q.  The  results  for  the  pulse  shapes 
are  shown  in  Figures  13  and  14,  relative  to  the  observed  results,  both 
unsealed  and  scaled.  It  appears  that  the  peak  values  of  the  pulses  are 
well  reproduced  but  the  peaks  are  smoothed  out  and  broadened  a  bit  so  that 
it  is  easily  possible  to  distinguish  the  two  cases. 

A  second  means  of  looking  at  a  Q  independent  of  amplitude  was 
carried  out  by  fitting  the  data  in  Figure  12  with  an  approximation  to  Q 
which  is  frequency  dependent  with  the  form 

Q  =  20f  (20) 
f  +  30 

although  this  has  no  corresponding  c  which  we  know  of  in  the  sense  of  a 
creep  function.  This  ad  hoc  Q  was  then  taken  with  a  constant  c  to  again 
propagate  the  initial  SALMON  pulse.  The  results  are  shown  in  Figures  15 
and  16  for  the  unsealed  and  scaled  cases  compared  to  the  original  data. 

In  this  case  of  a  frequency  dependent  Q,  the  peak  of  the  veloc¬ 
ity  pulses  follow  the  data  quite  well  and  there  is  a  strong  correspondence 
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Figure  15.  Estimates  of  pulses  generated  from  the  SALMON  record  at 
166  meters  using  the  fit  to  Q(«). 
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between  the  widths  and  sharpness  of  the  peaks  in  the  data  and  the  linear 
fit.  However,  It  Is  still  not  possible  to  provide  a  good  fit  to  the  ini¬ 
tial  portion  of  the  experimental  pulses  which,  as  Indicated  before,  are 
thought  to  be  due  to  shear  failure.  Still  the  remainder  of  the  pulses 
appear  to  be  consistent  with  a  Q  which  Is  Independent  of  amplitude. 

DISCUSSION 

We  have  taken  a  set  of  six  data  records  from  the  SALMON  event 
and  attempted  to  subject  them  to  such  an  analysis  that  any  Indication  of 
nonlinear  behavior  over  the  strain  regime  from  4  *  10"3  to  3  x  10“4  might 
be  detected.  Basically  we  find  that  the  SALMON  data,  taken  by  themselves, 
are  consistent  with  linear  but  Inelastic  behavior  which  can  be  described 
by  a  Q,  which  Is  Independent  of  amplitude  but  which  Is  at  least  a  mild 
function  of  frequency.  In  spite  of  the  fact  that  a  strain  variation  of 
more  that  an  order  of  magnitude  is  spanned  by  the  data,  the  closest  and 
most  distant  pairs  of  records  do  not  show  a  significant  difference  in 
effective  Q  values.  This  Q  approaches  a  value  of  about  20  at  frequencies 
above  100  Hertz  and  declines  to  zero  for  lower  frequencies. 

It  must  be  pointed  out  that  these  results  cannot  rule  out  the 
possibility  that  the  same  records  may  be  consistent  with  some  unspecified 
nonlinear  mechanism  which  simply  produces  effects  which,  locally,  are 
similar  to  linear  mechanisms.  In  fact,  this  Is  a  possibility  which  must 
be  taken  seriously  since  the  experiments  by  Tlttmann  suggest  that  for 
strains  of  less  than  10*6  a  Q  of  about  500  Is  appropriate.  Still  the 
experiments  are  under  conditions  sufficiently  different  from  explosions 
that  one  must  be  cautious  In  comparing  the  effects  under  somewhat  differ¬ 
ent  frequency  regimes  and  different  loading  conditions. 


SECTION  3 

ANALYTIC  APPROACH  TO  NONLINEAR  FRICTION 


In  discussing  the  departure  of  a  propagation  medium  from  purely 
elastic  behavior,  it  is  common  to  focus  on  the  quantity  Q,  the  absorption 
per  cycle.  If  one  wishes  to  calculate  the  detailed  propagation  of  a  seis¬ 
mic  wave,  which  is  what  we  shall  want  to  do  in  this  section,  one  needs 
constitutive  relations  which  can  be  combined  with  dynamical  laws  to  speci¬ 
fy  the  propagation.  For  a  linear  theory,  Q  and  the  phase  velocity  com¬ 
bine  together  to  specify  the  propagation;  they  contain  equivalent  infor¬ 
mation  to  constitutive  relations.  Furthermore,  Q  and  the  phase  velocity 
are  related  through  dispersion  relations.  For  a  nonlinear  theory,  both 
concepts,  Q  and  the  phase  velocity,  lose  precise  meaning  and  we  do  not 
know  how  to  translate  information  on  Q  into  (nonlinear)  constitutive  rela¬ 
tions.  While  not.  In  any  way  we  currently  understand,  providing  a  direct 
determination  of  the  constitutive  relations,  absorption  measurements  do 
provide  a  constraint  on  them  by  specifying  results  with  which  a  correct 
rheology  must  agree.  There  have  been  a  variety  of  efforts  to  relate  Q  to 
physical  properties  of  materials,  especially  internal  sliding  friction  due 
to  cracks.  Physically  based  models  of  attenuation  such  as  those  by 
Savage14,  for  thermoelastic  losses,  by  Mavko  and  Nur15  ,  for  viscous 
attenuation  by  partially  filled  cracks,  and  by  Walsh16  for  sliding  fric¬ 
tion  all  provide  a  means  of  estimating  the  loss  rate,  and  so  Q,  in  terms 
of  physical  parameters,  such  as  crack  length,  which  are  often  not  directly 
accessable  to  direct  measurement.  As  a  result  these  linear  models  provide 
Q  which  Is  only  loosely  tied  to  the  observables  of  the  problem  and  so  they 
do  not  provide  for  direct  verification  of  the  source  mechanism.  As  a 
result  these  models  are  not  much  more  satisfying  than  a  direct  phenomeno¬ 
logical  specification  for  Q. 


Many  experimental  studies6*  17» 18  of  attenuation  Indicate  that  It 
Is  generally  amplitude  dependent  with  an  onset  of  Increased  attenuation  at 
strains  of  about  10"6.  As  a  result  It  has  been  recognized  that  friction 
models  employing  mechanisms  which  have  a  nonlinear  basis  must  be  used  to 
account  for  the  data.  In  this  vein,  Mavko19  has  demonstrated  that  sliding 
friction  on  crack  surfaces  with  Irregularities  and  Including  partially 
closed  cracks  will  product  attenuation  which  Increases  with  amplitude. 
Using  a  Q  defined  defined  In  terms  of  relative  energy  loss  per  cycle, 
Mavko  finds  Q"1  proportional  to  the  peak  strain  over  the  cycle.  In  a 
similar  fashion,  Stewart  et  al.,20  have  proposed  a  physical  attenuation 
model  based  on  complex  crack  structure  and  friction  which,  like  Mavko, 
provides  a  Q-1  proportional  to  peak  strain  and  which  is  also  proportional 
to  the  crack  density  and  the  -4/3  power  of  the  confining  pressure. 

All  the  above  physical  models  for  friction  have  the  feature  that 
they  isolate  a  portion  of  the  material.  Impose  a  external  loading  and 
Integrate  the  energy  losses  over  a  cycle.  For  the  linear  case,  this  can 
be  translated  into  a  Q  function  which  can  be  Included  in  constitutive 
relations  which  will  apply  to  the  equations  of  motion  of  the  medium  and 
determine,  in  principle,  the  effects  of  attenuation  on  seismic  signals. 
However,  when  the  nonlinear  friction  models  Indicate  that  the  effective  Q 
is  a  function  of  peak  strain,  this  cannot  be  translated  directly  into  a 
corresponding  Q  since  Q  Is  basically  a  linear  concept.  As  a  result  there 
Is  currently  no  consistent  way  of  Including  the  nonlinear  attenuation  In 
calculations  of  seismic  propagation. 

In  this  section  we  wish  to  explore  a  possible  means  of  account¬ 
ing  for  nonlinear  attenuation,  at  least  In  a  phenomenological  way  which  Is 
consistent  with  the  data  and  which  Is  not  confined  to  effective  Q  Inte¬ 
grated  over  a  cycle  assumptions.  This  will  be  done  by  consideration  of 
constitutive  relations  which  are  explicitly  nonlinear  and  which  allow 
approximate  analytic  solutions  so  that  their  consequences  can  be  readily 
explored. 


We  shall  specify  constitutive  relations  which  provides  for  pro¬ 
pagation  which  satisfies  Items  1,  2,  and  3  given  In  the  introduction. 

1)  The  absorption  of  energy  Increases  with  ampli¬ 

tude;  one  relation  which  has  been  proposed  to 
express  this  effect  is: 


l*  V  ' 


where  e  is  the  peak  strain;  the  quantity  Q  will 
be  discussed  below. 

The  propagation  satisfies  cube-root  scaling;  that 
Is,  the  propagation  Is  invariant  under  the  scale 
transformations: 

t  ♦  y1  /3  ts 
x  +  y1  /3  xj 
u  -*•  y1  /3  us 


r-V-V-'. 


where  t$,  x$  and  u$  are  respectively  the  scaled 
time,  scaled  distance  and  scaled  displacement 
while  y  Is  the  event  yield. 


3)  The  propagation,  although  presumably  nonlinear 
(property  1),  nevertheless  provides  a  waveform 
which  changes  only  slowly  with  distance. 

In  addition  to  providing  for  propagation  with  these  properties 
the  relations  we  shall  use  are  Galilean  Invariant  and  satisfy  causality — 
In  fact  the  relations  are  local  in  both  space  and  time. 


-‘'4 
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INTRODUCING  FRICTION 


We  start  from  the  elastic  equation  in  one  cartesian  dimension: 

p  aiy.  -  u  lii  (21) 

at?  3X2 

where  u  is  the  displacement  and  u  is  an  elastic  modulus.  We  can  introduce 
friction  into  the  system  by  allowing  the  modulus  to  change  with  the  strain 
and  rate  of  change  of  strain  in  such  a  way  that  the  increment  in  force 
always  opposes  changes  in  the  state  of  the  system: 

u  +  u(  1  +  o  sgn (<LH  .))  (22) 

3  X  3  X3 1 

A  form  similar  to  (22)  Is  given  by  Knopoff  and  McDonald21;  we  prefer  the 
form  (22)  given  It  is  Galilean  invariant  whereas  that  of  Reference  21  is 
not;  for  propagating  disturbances  -  those  for  which  the  x  and  t  deriva¬ 
tives  are  related  by  a  constant  -  the  two  forms  are  equivalent.  In  Refer¬ 
ence  21  it  is  argued  by  analogy  with  single  oscillator  models  that  (22) 
will  lead  to  a  frequency  independent  Q.  By  this  same  argument  (22)  will 
lead  to  an  amplitude  Independent  Q. 

To  satisfy  our  item  1  -  that  absorption  Increases  with  amplitude 
-  we  make  the  "coefficient  of  friction"  a  depend  on  the  strain.  One  pos¬ 
sibility  is: 


a 


(23) 


In  this  case  a  Increases  with  material  distortion  whether  the  material  is 
under  tension  or  under  compression.  Such  need  not  be  the  case  and  we 
might  equally  well  have: 


Which,  if  any,  of  these  relations  is  appropriate  for  materials  in  the 
earth  is  an  experimental  question.  Lacking  such  guidance  we  shall 
consider  relation  (23);  the  qualitative  conclusions  reached  here  would 
apply  to  the  other  relations  as  well  although  the  details  of  the 
propagation  would  be  different. 


We  have  arrived  at  the  following  dynamical  equation: 


It  is  easy  to  see  that  this  equation  is  yield1/3  scale  invariant,  that  is, 
it  Is  form  invariant  under  the  scale  transformations  of  item  2. 

APPROXIMATE  SOLUTIONS 

In  this  report  we  are  primarily  concerned  with  amplitude 
dependent  effects.  The  number  y  in  relation  (23)  determines  the  amplitude 
independent  value  the  absorption  will  approach  in  the  limit  of  small 
amplitudes.  To  simplify  future  expressions  and  clarify  the  amplitude 
dependent  effects  predicted  by  equation  (23)  we  shall  set  y  equal  to 
zero.  Furthermore  we  stated.  In  Item  3,  that  we  wished  to  have 
propagation  of  an  approximately  unchanging  waveform.  Thus  we  shall  seek 
solutions  of  the  form: 


u  =  F(x)  G(t  -  X) 
c 


where  c  is  a  constant  and  we  shall  normalize  the  functions  F  and  G  by 
specifying  Max  |g|  *  1.  Substitution  of  (27)  into  (26)  does  not,  without 
further  assumptions,  yield  a  simple  expression.  If,  however,  we  further 


assume  that: 


Fu 
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where  Xg  is  the  length  scale  of  the  waveform  G,  that  is,  we  assume  that 
the  fractional  change  in  the  amplitude  function  F  when  the  wave  propagates 
through  its  length  is  small,  then  equation  (26)  reduces  to: 


F_L  =  -  ofi_  |g" 

F2  2c2  1 


where  we  have  chosen  c  *  /p/p  ;  the  wave  travels  at  the  linear  sound 
speed.  The  left  hand  side  of  (29)  depends  only  on  x  whereas  the  right 
hand  side  depends  only  on  the  retarded  time  t  (t  it-  *.)•  Thus  each  side 
must  be  equal  to  a  constant  and  G  is  a  quadratic.  Thus  in  somewhat  of  an 
analogy  to  soliton  theory  there  is  only  one  class  of  waveforms 
which  will  propagate  according  to  (27). 

The  general  solution  for  G  is  given  by: 

G  *  4  t (T-t i  0  <  t  <  T  (30) 

l2 

where  T  is  a  constant  -  the  duration  of  the  wave.  From  this  expression  we 


G" |  =  _S 
1  T2 
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Combining  this  expression  with  (29)  we  find  the  general  solution  for  F: 


F  =  I  1 
4a0  a  +  x 


where  a  is  an  arbitrary  constant.  The  displacement  is  thus  given  by: 

„  =  c?T  1  t(T-t)  (33) 

o$  a  +  x  T2 

For  a  propagating  linear  wave  of  frequency  u>  a  common  definition 
of  Q  is  in  terms  of  the  log-derivative  of  the  amplitude;  if  A  is  the 
amplitude  we  have: 

d_  J.pA  =  -  Q-i  (34) 

dx  2c 

Since  F  is  the  amplitude  of  our  nonlinear  wave,  the  natural  definition  of 
Q  is 

QJ  =  -  2c  £!  =  &£_  p  (35) 

til  F  u)  cT2 

Thus  we  see  that  for  the  current  model  Q"1  is  a  linear  function  of  the 
amplitude  as  Minster  and  Archambeau  suggested  might  be  characteristic  of 
the  SALMON-COWBOY  data.  The  characteristic  frequency  of  our  wave  will  be 


0)  = 

T 


which  gives: 


QJ  =  SaBm  F 


This  relation  demonstrates  our  earlier  statement  that  for  nonlinear 
propagation  cube-root  scaling  does  not  imply  the  frequency  independence  of 
Q:  equation  (26)  is  scale  invariant  but  Q  depends  on  id. 


given  by: 


From  equation  (33)  we  find  that  the  peak  velocity  in  the  wave  is 


v  =  cl  l 
P  as  a  +  x 


This  relation,  and  all  our  work  so  far  applies  only  to  propagation  in  one 
cartesian  dimension.  Propagation  in  one  spherical  dimension  will  be  more 
complicated  and  we  have  not  yet  carried  out  a  proper  analysis  of  that 
case.  If  we  concern  ourselves  with  radii  which  are  not  too  small,  so  that 
the  afterflow  term  is  small,  we  expect  that  the  peak  velocity  will  behave 
like: 

v  -  1  (39) 

P  08  r(a  +  r) 

Recalling  the  peak  velocity  data  from  SALMON,  we  had  vp  =  14  m/s  at  166 
meters  and  vp  =  1  m/s  at  660  meters.  Using  these  two  data  points  to 
determine  the  constants  in  (39)  we  find  *  4.6  x  1$  m3  /s  and  a  =  30 
meters.  For  other  stations  in  the  SALMON  data  we  get 

Table  2.  Comparison  of  fit  with  SALMON  data. 
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It  appears  from  the  table  that  we  get  a  pretty  good  fit  to  the 
SALMON  peak  velocity  data  using  the  nonlinear  friction  model.  As  we  saw 
in  the  previous  section  however,  the  SALMON  shock  propagated  not  according 
to  the  nonlinear  friction  model  but  rather  as  a  wave  in  an  analastic 
solid.  We  do  not  know  what  to  make  of  the  apparently  good  fit  shown  in 
the  table;  probably  one  should  view  it  as  an  example  of  the  danger  inher¬ 
ent  in  drawing  too  firm  a  conclusion  about  the  validity  of  propagation 
models  using  only  peak  velocity  data. 

DISCUSSION 

We  have  seen  that  the  constitutive  relation  given  in  equation 
(23)  provides  for  propagation  having  the  properties  set  out  at  the  begin¬ 
ning  of  this  chapter.  While  we  know  from  the  detailed  analysis  of  the 
SALMON  data  given  in  the  previous  chapter  that  the  SALMON  shock  did  not 
possess  all  these  properties  -  in  particular  item  1  -  it  may  be  that  the 
COWBOY  data,  taken  at  lower  strains,  will  exhibit  the  properties  and  that 
the  nonlinear  friction  model  will  be  useful  in  modeling  these  results.  We 
will  not  know  whether  or  not  this  is  the  case  until  we  complete  a  detailed 
study  of  the  COWBOY  data  in  the  same  way  we  have  done  for  SALMON. 


SECTION  4 
CONCLUSIONS 


From  the  studies  presented  in  this  report  we  draw  the  following 
principal  conclusions: 

1.  The  seismic  wave  from  event  SALMON  propagated,  in  the  range 
of  strains  from  10“3  to  10-4  ,  very  much  like  a  wave  in  an 
anelastic  solid. 

2.  At  some  point  in  the  propagation  the  attenuation  must 
become  less  in  order  to  join  smoothly  to  small  strain 
measurements. 

3.  By  applying  the  procedures  used  in  the  present  report  for 
SALMON  to  the  COWBOY  data  we  should  be  able  to  study  the 
nature  of  the  transition  in  some  detail. 

4.  Nonlinear  friction  can  provide  for  propagation  with  proper¬ 
ties  much  like  those  thought  to  be  characteristic  of  the 
body  of  data  for  intermediate  strains  (though  not  for 
SALMON). 

The  arguments  for  item  2  are  that  Tlttmann's  laboratory  data  taken  at 

strains  of  about  10"6  show  much  smaller  absorption  and  that  propagation 
through  a  medium  as  highly  absorbing  as  we  measure  for  SALMON,  would  not 
reach  teleselsmic  distances;  while  it  Is  true  that  teleseismic  paths 
entirely  through  salt  are  not  observed  It  Is  not  generally  thought  that 
natural  halite  formations  have  low-amplitude  absorption  so  much  greater 
than  other  seismic  media. 
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ANALYTIC  APPROACHES  TO  LINEAR  AND  NONLINEAR  ATTENUATION* 


SUMARY  A  review  of  near-field  data  from  explosions  indicates  that 

nonlinear  behavior  extends  well  beyond  a  few  hundred 
meters/kt1 /3 .  Frequency  dependent  contributions  to  Q, 
perhaps  from  scattering,  can  lead  to  underestimates  of 
magnitudes  using  spectral  ratio  methods. 

INTRODUCTION 

Until  recently  it  has  been  assumed  that  the  behavior  of  rock 
beyond  the  radius  of  cracking  (>  100  m/kt1 /3 )  due  to  a  explosive  source 
shows  a  linear  but  perhaps  inelastic  character.  Beyond  -300  m/kt1/3,  the 
"elastic  radius,"  seismic  propagation  has  been  taken  as  nearly  elastic 
with  mild  attenuation  due  to  anelasticity  which  is  generally  described  in 
terms  of  Q  or  Q“l  ,  the  "internal  friction."  As  a  result  a  major  topic  for 
study  for  seismic  detection  is  that  of  the  behavior  of  path  integrated  Q-1 
as  a  function  of  path,  frequency  and  possibly  amplitude.  The  approach  to 
the  study  of  QJ  has  been  largely  phenomenological  and  experimental  with 
only  moderate  consideration  of  physical  mechanisms  or  use  of  analytic 
techniques. 

In  this  paper  we  shall  review  the  data  from  explosions  at  ranges 
from  about  10  m/kt1/3  to  3000  m/kt1/3  in  order  to  establish  the  extent  to 


♦This  paper  was  submitted  to  the  DARPA-sponsored  meeting  on  BASIC  RESEARCH 
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which  the  attenuation  must  be  described  as  nonlinear  and  the  extent  to 
which  Q"1  Is  dependent  upon  frequency.  It  will  be  found  that  these  near¬ 
field  data  suggest  a  significant  nonlinearity  for  which  an  empirical  Q"1 
is  consistent  with  being  linearly  dependent  upon  amplitude  at  large  strain 
but  which  becomes  constant  at  small  strains. 

A.  NEAR-FIELD  DATA  AND  NONLINEAR  ATTENUATION 

Near-field  data  from  explosive  sources  in  salt  provides  the  best 
available  information  on  the  attenuation  of  pulses  in  the  moderate  strain 
regime.  These  data  are  from:  Larson's  work  in  pressed  salt  with  scaled 
ranges  from  10  m/kt1 /3  to  200  m/kt1/3  {mean  frequency  -  id*  Hz);  Salmon 
nuclear  explosion  in  dome  salt  with  ranges  from  100  m/kt1/3  to  400  m/kt1/3 
(mean  frequency  -5  Hz);  and  Cowboy  series  tamped  chemical  explosions  in 
dome  salt  with  ranges  from  200  m/kt1/3  to  3000  m/kt1/3  (mean  frequency  -10 
Hz  or  greater).  A  plot  of  peak  velocity  versus  scaled  range  is  shown  in 
Figure  A-l.  This  illustrates  the  well  known  result  that  simple  cube  root 
scaling  applies  with  remarkable  accuracy  to  data  from  events  with  a  yield 
variation  of  ten  orders  of  magnitude  (and  thus  a  range  of  frequencies 
spanning  more  than  three  orders  of  magnitude).  This  result  is  all  the 
more  Interesting  because  It  takes  place  In  a  regime  where  the  peak  values 
of  both  velocity  and  displacement  fall  off  much  more  rapidly  than  1/r. 
That  is,  for  this  spherical  geometry,  the  behavior  cannot  be  simple 
elastic  since  the  decay  exceeds  that  for  geometrical  divergence  alone. 
The  residual  attenuation  may  result  from  linear  but  inelastic  (l.e.,  an- 
elastic)  behavior  or  from  nonlinear  attenuation.  In  either  case,  the 
exi stance  of  simple  scaling  demonstrates  that  the  effective  QJ  must  be 
nearly  Independent  of  frequency,  at  least  for  this  frequency  and  strain 
regime.  The  peak  strains  (-  peak  velocity/propagation  speed)  from  Figure 
1  range  from  about  10"1  to  10-5  . 
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Figure  A-l.  Peak  particle  velocities  from  explosions  In  salt, 
from  Larson  (1982). 


In  order  to  address  the  question  of  linear  or  nonlinear  attenua¬ 
tion,  we  now  consider  the  implications  for  such  possibilities.  If  a 
linear  viscoelastic  mechanism  Is  to  account  for  these  data,  the  effect  can 
be  expressed  in  terms  of  a  Q  which  is,  at  least  approximately,  independent 
of  frequency  over  the  frequency  range  of  the  experiments.  This  leaves 
some  latitude  in  selection  of  an  attenuation  model  but  since  the  behavior 
In  other  frequency  regimes  has  little  effect,  a  reasonable  attempt  Is  to 
simply  take  Q  to  be  completely  Independent  of  frequency.  Such  a  model  has 
been  constructed  by  Kjartansson  (1980)  who  provides  rather  simple  expres¬ 
sions  for  attenuation  and  dispersion  which  can  be  written  as  a  complex 
phase  velocity 


:in_  sgnU ) 


C(dJ  ) 


where  ^  is  a  modulus,  p  the  density,  W()  an  arbitrary  frequency  scale  and 
y  is  an  arbitrary  measure  of  the  attenuation  rate.  In  particular,  y  « 
Q"1  /tt  if  Q  >  1/jr  -  For  a  spherical  pulse  with  a  velocity  spectrum  at  radi¬ 
us  'a'  of  V(a rfo )  the  velocity  pulse  at  any  time  is  given  by 
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where 


is  the  phase  velocity 


and 


o  *  (o  c-1  tan(»y  sgn{«o)/2)  is  the  spatial  attenuation 
constant. 


A  velocity  pulse  from  Salmon  at  a  scaled  range  of  230  m/kt1 /3  is 
shown  by  Larson  (1982)  and  this  is  used  as  the  initial  point  for  a  study 
of  effects  of  anelastic  attenuation  for  constant  Q.  A  sequence  of  pulses 
at  increasing  range  is  shown  in  Figure  A-2  Including  frequency  dependent 
attenuation  (a  is  proportional  to  u>  for  constant  Q)  and  the  attendant  dis¬ 
persion.  Attenuation  and  lengthening  of  the  pulse  are  evident  for  this 
case  using  Q*10.  The  corresponding  decay  of  peak  velocity  as  compared  to 
the  Salmon  and  Cowboy  data  Is  shown  In  Figure  A-3.  The  value  of  Q  was 
chosen  to  match  approximately  the  slope  of  the  Salmon  data  (vmflX  '  r_1  *88 ) 
at  the  Initial  range  and  It  Is  apparent  that  at  larger  ranges  the  decay  is 
somewhat  more  rapid  for  constant  Q.  (U  is  easy  to  show  that  for  suffi¬ 
ciently  large  ranges  the  peak  velocity  will  fall  off  like  r-2  if  the  spec¬ 
trum  is  flat  below  the  corner  frequency.  Smaller  Q  (greater  attenuation) 
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Figure  A-2.  Waveform  of  SALMON  velocity  pulse  at  three  ranges  • 
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results  in  a  more  rapid  onset  of  the  r-2  behavior.)  From  this  it  seems 
that,  in  the  absence  of  other  constraints,  a  rather  small  Q  -  10  from  an 
anelastic  model  can  give  a  rough  fit  to  the  Salmon  and  Cowboy  data.  How¬ 
ever  recent  work  at  Rockwell  (Tittmann,  1983)  indicates  that  for  strains 
from  10"8  to  10"6 ,  the  Q  from  dome  salt  is  approximately  500  as  measured 
at  -  500  Hz.  Tittmann  further  notes  that  at  strains  greater  than  10-6  Q"1 
increases  significantly  at  least  up  to  strains  of  10-5  which  is  the  maxi¬ 
mum  available  in  the  experiment.  In  short  the  experiments  suggest  that 
attenuation  increases  with  amplitude  (i -e. ,  it  is  nonlinear)  for  strains 
in  excess  of  10”6  •  The  magnitude  of  the  measured  Q  at  small  strains  is 
much  different  from  that  required  to  fit  the  explosion  data  suggesting 
that  a  linear  mechanism  cannot  account  for  these  data.  Furthermore  the 
Rockwell  data  for  nonlinear  behavior  indicate  an  onset  at  - 10"6  suggesting 
that,  since  even  the  smallest  strain  from  Cowboy  is  - 10-5  ,  nonlinear  be¬ 
havior  is  a  vital  element  in  accounting  for  the  decay  of  the  explosively 
generated  pulses  in  the  near-field.  We  intend  to  apply  the  spectral  ratio 
method  to  Salmon  and  Cowboy  data  to  verify  nonlinearity  as  well  as  to 
explore  any  frequency  dependence  of  Q. 

Microscopic  models  of  frictional  attenuation  due  to  sliding  of 
crack  surfaces  such  as  Mavko  (1980)  and  Stewart  et  al .  (1982)  generate  a 
contribution  to  Q_1  which  is  proportional  to  the  peak  strain  during  a 
cycle  but  which  is  also  independent  of  frequency.  As  a  result  it  has  been 
suggested  that  an  approximate  Q  can  be  written  as 

Q-1  -  Off1  +  (e/e„)  Q;1 

eo 

where  is  the  "intrinsic"  Q  and  the  second  term  is  the  frictional  Q 
which  is  linear  in  strain  (and  normalized  at  e0  ).  Such  a  form  has  been 
considered  by  Minster  and  Archambeau  (1983)  in  a  fashion  similar  to  that 
which  follows  here.  The  above  form  for  a  nonlinear  Q  ignores  the  fact 
that  Q  is  Inherently  a  linear  concept.  Here  it  is  to  be  interpreted  to 
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mean  that  the  spectral  content  will  decay  with  range,  r,  according  to  a 
factor 


exp(  -  ^JL) 
^  ZcQJ 


where  the  strain,  e,  which  appears  is  taken  as  the  peak  strain  during  the 
pulse.  Using  such  a  model  for  nonlinear  behavior,  we  have  found  the 
attenuation  of  the  Salmon  velocity  pulse  by  a  closely  spaced  series  of 
Fourier  synthesis  advancements  in  which  Q-1  iS  taken  to  be  constant  in 
each  small  step  but  altered  (in  fact  decreased)  according  to  the  amplitude 
between  steps.  The  limiting  result  for  sufficiently  small  steps  is  shown 
in  Figure  3.  In  the  example  shown  t^j  =103  and  Q  at  the  smallest  range  (230 
m/kt1  /3  )  was  5.  These  variable  Q  parameters  have  been  chosen  both  to  give 
a  reasonable  fit  to  the  Salmon  data  (which  are  always  a  bit  above  the  Cow¬ 
boy  data)  and  to  provide  an  onset  of  nonlinear  behavior  which  is  consis¬ 
tent  with  the  Rockwell  experiments.  However  the  result  is  rather  insensi¬ 
tive  to  the  value  of  because  the  smallest  strain  in  the  explosion  data 
is  wel 1  above  10-6  . 


It  appears  that  a  model  of  nonlinear  Q-1  with  a  dominant  term 
proportional  to  the  strain  gives  a  decay  of  the  pulse  from  explosions  in 
salt  which  is  consistent  with  observation  as  well  as  consistent  with  small 
strain  behavior  of  Q.  However  the  data  are  not  adequate  to  fix  the  expon¬ 
ent  of  the  strain  in  Q”1  .  For  example  Q"1  «  e1  or  e2  can  also  give  a 
fair  fit  to  the  data  although  the  transition  to  linear  behavior  is  moved 
somewhat  from  the  10-6  value.  The  strongest  statement  we  can  make  is  that 
the  assumed  form  of  Q  is  consistent  with  the  explosion  3nd  laboratory  data 
but  the  particular  form  is  not  tightly  constrained.  Still  a  nonlinear 
mechanism  appears  to  be  necessary  to  account  for  the  two  types  of  data. 


B. 


IMPLICATIONS  OF  Q  BEHAVIOR  FOR  TELESEISMIC  SIGNALS 


Customary  calculations  of  teleseisimic  propagation  of  signals 
from  explosions  and  earthquakes  often  treat  a  near-source  region  as  giving 
the  initial  conditions  for  a  subsequent  linear  propagation.  In  the  case 
of  explosions  this  "elastic  radius"  is  often  taken  as  -300  m/kt1/  3  and  it 
is  assumed  that  all  nonlinear  behavior  occurs  inside  this  radius.  It  is 

now  generally  though  that  nonlinear  effects  extend  out  to  larger  radii, 
perhaps  out  to -10*  m/kt1/3.  The  behavior  in  this  intermediate  strain 
regime  is  a  significant  source  of  uncertainty  in  the  ultimate  problem  of 
relating  teleseismic  signals  to  source  characteristics. 

From  the  point  of  view  of  linear  calculations  of  seismic  sig¬ 
nals,  nonlinear  mechanisms  provide  anomolous  attenuation  which  has  the 
potential  for  altering  the  relation  between  teleseismic  signals  and  the 
near-field  pulse.  However  mere  existance  of  nonlinear  attenuation  does 
not  necessarily  negate  spectral  ratio  methods  unless  the  nonlinear  aspect 

has  a  different  frequency  dependence  that  the  assumed  Q.  This  is  due  to 

the  fact  that  spectral  ratio  methods  depend  upon  finding  a  path-integrated 
q-1  ,  which  is  generally  expressed  at  t*,  which  may  or  may  not  be  a  func¬ 
tion  of  frequency.  If,  for  example,  t*  is  taken  as  independent  of  fre¬ 
quency  and  found  by  matching  the  high-frequency  slope  of  the  observed 
spectrum  to  an  assumed  source  spectrum,  the  existance  of  a  frequency 

independent  but  nonlinear  contribution  to  t*  would  not  alter  the  result 

except  to  raise  the  value  of  t*  which  is  proportional  to  the  mean  value  of 
Q-1  along  the  path.  On  the  other  hand,  if  the  nonlinear  contribution  to 
t*  has  a  frequency  dependence  different  from  the  assumed  one,  use  of  a 
suitable  slope  connecting  t*  to  find  a  magnitude  at  the  source  may  be  in 
error.  This  question  has  been  raised  by  Der  (1983)  in  a  slightly  differ¬ 
ent  context.  This  arises  from  use  of  spectral  shape  changes  with  propaga¬ 
tion  distance  to  deduce  information  about  spectral  amplitude  changes.  In 
particular  It  is  noted  that  If  the  true  t*,  which  describes  the  amplitude 
change,  is  related  to  the  spectral  slope  apparent  t*  by 


t*  =  t*  +  (i)  3  t* 
do) 

If  t*  is  independent  of  u  the  two  are  identical  but  if  t*  decreases  with 
frequency  t*<t*.  According  to  the  work  of  Lundquist  and  Samowitz  (1983), 
such  is  in  fact  the  case  although  this  remains  a  question  for  discussion. 
For  example,  Der  and  Lundquist,  working  with  exactly  the  same  data  set, 
have  found  significantly  different  values  for  t*.  Using  an  absorption 
band  model  for  Q  (for  which  Q  «  u  at  large  u)  Lundquist  has  fit  the  data 
using  as  variables  the  low-frequency  t*  and  a  period  below  which  r*  rolls 
off.  Der,  on  the  other  hand,  has  fit  the  data  assuming  t*  is  constant  in 
various  frequency  bands  and  then  inverted  the  above  relation  to  find  t*. 
At  2Hz  the  difference  between  the  fits  at  distant  stations  appear  to  be 
At*  “.1  to  .2  sec  so  this  translates  into  a  source  intensity  ratio  of  ”10. 
That  is,  attenuation  will  be  underestimated  if  t*  is  fit  as  constant  and 
this  will  cause  the  magnitude  of  the  source  to  be  underestimated. 

It  should  be  noted  that  if  QJ  has  a  term  which  is  proportional 
to  u”1  this  will  serve  to  provide  an  overall  attenuation  of  the  spectral 
amplitude  but  it  will  not  alter  the  spectral  shape  so  that  the  contribu¬ 
tion  to  t*  will  be  zero.  Thus  it  is  very  important  to  include  any  such 
effects  which,  for  example,  occur  at  high  frequency  in  the  absorption  band 
model.  Observations  and  analysis  by  Aki  and  Chouet  (1975)  of  coda  attenu¬ 
ation  indicate  that  at  least  some  portions  of  the  coda  result  from  single 
scatterings  of  the  main  pulse  by  inhomogeneities  in  the  lithosphere. 
Earthquake  coda  have  been  analyzed  to  find  an  effective  Q  which  is  fre¬ 
quency  dependent  giving  Q  increasing  with  frequency  from  1  to  25  Hz. 
Dainty  (1983)  has  looked  at  other  earthquakes  finding  similar  effects 
showing  Q,  apparently  due  to  scattering,  which  is  approximately  propor¬ 
tional  to  frequency  in  the  3  to  50  Hz  regime. 


Dainty,  in  a  theoretical  analysis,  has  found  the  contribution  of 
scattering  to  Q  by  considering  the  Born  approximation  in  the  presence  of 
randomly  distributed  inhomogeneities  in  propagation  speed  in  the  medium. 
It  is  found  that  if  the  autocorrelation  function  of  the  inhomogeneities 
has  a  nonzero  first  derivative  at  zero  lag,  the  contribution  to  Q  is  pro¬ 
portional  to  frequency.  It  is  suggested  that  discontinuities  at  the  scat¬ 
tering  structure  edge  lead  to  this  result.  Basically  the  Born  scattering 
Q  will  be  proportional  to  u  for  wavelengths  small  compared  to  the  overall 
dimensions  of  the  structures  but  large  compared  to  the  thickness  of  the 
edges. 


Independent  of  the  exact  nature  of  the  physical  mechanism  which 
leads  to  the  lithospheric  contribution  to  Q,  the  fact  that  the  observed 
character  is  Q«  u  Is  of  interest  for  application  of  spectral  ratio  meth¬ 
ods  to  waves  which  pass  through  the  lithosphere.  Both  Aki  and  Dainty 
indicate  that  the  Q  found  from  coda  takes  a  value  of  “200  at  frequencies 
of  a  few  Hertz.  As  was  previously  noted  a  Q_I  *  w"1  will  provide  an  over¬ 
all  attenuation  of  the  amplitude  of  the  pulse  spectrum  but  it  will  not 
alter  the  shape.  An  underestimate  of  source  magnitude  will  result  from 
propagation  through  the  lithosphere  and  the  degree  can  be  found,  given  a 
path  length.  For  example  for  a  path  of  100  km  an  additional  Q-1  of  5  x 
10-3  at  2  Hz  will  lead  to  an  apparently  decreased  source  magnitude  at  all 
frequencies  of 


exp(  ID  -  5 
cQ 

if  the  source  of  attenuation  is  not  explicitly  taken  into  account. 

C.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  preliminary  study  of  the  attenuation  observed  for  explosively 
generated  near-field  pulses  in  salt  indicates  that  the  behavior  is  dis- 


tinctly  nonlinear  for  the  strains  for  10-1  to  10"5 .  The  observed  behavior 
is  consistent  with  a  Q"1  which  is  the  sum  of  a  constant  (Q  at  small 
strains  are  found  in  Rockwell  experiments)  and  a  term  linear  in  the  strain 
(which  is  as  found  in  frictional  crack  models).  It  is  recommended  that 
dynamical  models  be  formulated  which  conform  to  these  near-field  data  for 
inclusion  in  calculations  of  waveforms  out  to  a  proper  elastic  radius. 


A  frequency  dependent  Q  will  alter  the  relation  between  spectral 
amplitude  and  spectral  shape  leading  possibly  to  substantial  underesti¬ 
mates  of  source  magnitudes  using  spectral  ratio  methods.  This  is  espe¬ 
cially  true  for  lithospheric  scattering  which  gives  a  contribution  to  Q 
which  may  not  alter  the  spectral  shape  at  all.  It  is  important  to  estab¬ 
lish  the  extent  and  character  of  scattering  contributions  to  Q  in  order  to 
refine  magnitude  estimation  methods. 
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APPENDIX  B 
SALMON  DATA 


SALMON  was  a  5.3  kT  nuclear  explosion  which  took  place  at  a 
depth  of  800  meters  in  a  salt  dome  in  Mississippi  on  October  22,  1964.  It 
was  part  of  a  decoupling  experiment  in  which  explosion  Sterling  took  place 
In  the  SALMON  cavity  at  a  later  time.  The  SALMON  explosion  was  one  of  the 
best  instrumented  cases  ever  carried  out  due  to  the  unique  character  of 
the  experiment.  As  a  result,  the  data  from  SALMON  have  proved  to  be  the 
best  available  for  nuclear  explosions  in  terms  of  the  intermediate  range 
ground  motion  recordings.  Furthermore,  by  virtue  of  the  rather  uniform 
character  of  the  salt  dome,  the  data  are  internally  quite  consistent. 

The  data  were  originally  recorded  by  analog  methods  by  Sandia 
and  they  were  later  converted  to  digital  form  at  that  laboratory.  The 
Intermediate  range  data  from  166  meters  out  to  660  meters  were  analyzed 
in  great  detail  by  Perret1  as  well  as  by  Rogers2  and  Patterson3  shortly 
after  the  event.  In  more  recent  times  some  aspects  of  the  data  have  been 
reexamined  by  Murphy4,  Perl5,  and  TrulioP,  and  informally  by  others.  The 
taped  data  records  reported  here  were  supplied  to  us  by  Neil  Perl  of 
Applied  Theory,  Inc.,  who  originally  received  them  from  Sandia. 

Figure  8-1  indicates  the  arrangement  of  sensors  for  the  SALMON 
event.  Most  of  the  Instrument  stations  had  several  sensors  so  there  is 
some  redundancy  in  the  records  as  well  as  Information  concerning,  often, 
all  three  axes  of  motion.  The  records  supplied  to  us  are  of  two  general 
types,  acceleration  (A)  and  velocity  (U).  The  orientations  are  either 
horizontal  radial  (R),  vertical  (V)  or  horizontal  transverse  (T).  The 


records  available  are  given  in  Table  B-l  with  the  designation  EI-D-AO 
where  #  indicates  the  bore  hole  number,  D  the  depth  in  hundreds  of  feet, 
A  acceleration  or  velocity  and  0  the  orientation. 


Table  B-l.  SALMON  data  records. 


E5-S-AR 

E6-27-UR 

El 1-34 -AR 

E14C-27-AR 

-S-AT 

-27-UT 

-34 -AV 

-27-AV 

-S-AV 

-27-UV 

-34 -UR 

-27-UV 

-27-AT 

Ell-S-AR 

-34 -UV 

-32-AR 

-27-AV 

-S-AT 

E14-20-AR 

-32-AV 

-27-UT 

-S-AV 

-20-AV 

-32-UR 

-27-UV 

-S-UR 

-20-UR* 

-32-UV 

E6-S-AR 

-S-UV 

-34 -AR 

-36-AR 

-S-AT 

-20-AR 

E14C-S-AR 

-36-AV 

-S-UR 

-2Q-UR* 

-S-AV 

-36-UR 

-S-UV 

-2Q-UV 

-S-UR 

-36-UV 

-20-AR 

-27-AV 

-22-AR 

-39-AR 

-20-AV 

-27 -UR 

-22-AV 

-39-UR 

-27-AR 

-27-UT 

-22-UR 

-39-UV 

-27-AV 

-27-UV 

-22-UV 

*  Instrument  Mis-oriented 


The  records  of  interest  in  our  study  are  the  subsurface  cases 
excluding  transverse  motion.  The  transverse  motion  measurements  would  be 
zero  for  perfect  spherical  symmetry  and  for  perfect  orientation.  (Of 

course,  they  are  not  exactly  zero.)  We  have  looked  only  at  the  AR,  AV, 

UR,  and  UV  cases  excluding  those  at  the  surface.  The  records  consist  of 
approximately  27,700  pairs  of  samples  at  intervals  of  0.2  milliseconds. 
Each  pair  consists  of  a  time  in  units  of  seconds  and  an  acceleration 
(units  of  g,  the  acceleration  of  gravity)  or  a  velocity  (units  of  feet  per 

second).  The  time  ranges  from  about  -0.544  seconds  to  5.0  seconds  at 

Integer  multiples  of  200  microseconds. 


E527UV  SALMON 


E620AR  SALMON 


Figure  B-2.  Raw  SALMON  velocity  records 
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E627UR  SALMON 


E627UV  SALMON 


NMt(5EC) 

T)ME(SEC) 

E1120AR  SALMON 

E1120UR  SALMON 

.2  1.4  1.4  1.4  '  B.O  0.2  0.4 

0.4  0.4  1.0 

Figure  B-2.  (Continued). 
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Figure  B-2.  (Continued). 
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Figure  B-2.  (Continued) 
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